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Abstract — A two-dimensional grid with dots is called a config- 
uration with distinct differences if any two lines which connect 
two dots are distinct either in their length or in their slope. 
These configurations are known to have many applications 
such as radar, sonar, physical alignment, and time-position 
synchronization. Rather than restricting dots to lie in a square 
or rectangle, as previously studied, we restrict the maximum 
distance between dots of the configuration; the motivation for this 
is a new application of such configurations to key distribution 
in wireless sensor networks. We consider configurations in the 
hexagonal grid as well as in the traditional square grid, with 
distances measured both in the Euclidean metric, and in the 
Manhattan or hexagonal metrics. 

We note that these configurations are confined inside maximal 
anticodes in the corresponding grid. We classify maximal anti- 
codes for each diameter in each grid. We present upper bounds on 
the number of dots in a pattern with distinct differences contained 
in these maximal anticodes. Our bounds settle (in the negative) 
a question of Golomb and Taylor on the existence of honeycomb 
arrays of arbitrarily large size. We present constructions and 
lower bounds on the number of dots in configurations with 
distinct differences contained in various two-dimensional shapes 
(such as anticodes) by considering periodic configurations with 
distinct differences in the square grid. 

Index Terms — Anticodes, Costas Arrays, Distinct-Difference 
Configurations, Golomb Rectangles, Honeycomb Arrays 

I. Introduction 

A Golomb ruler (or ruler for short) of order m (also known 
as a Sidon set) is a set S of integers with jS"! = m having 
the property that all differences a — b (for a, 6 G S, with a 7^ 6) 
are distinct. They were first used by Babcock, in connection 
with radio interference [1]. The length of a Golomb ruler S 
is the largest difference between any two elements of S. It 
is easy to show that a ruler of order m has length at least 
(™); a ruler meeting this bound is called perfect. Golomb 
has shown that no perfect ruler exists with order greater than 
four [2]. The problem of finding the shortest possible length 
of a Golomb ruler of a given order has been widely studied; 
no general solution is known, but optimal rulers have been 
determined for orders less than 24 (see [3] for details). The 
elements of a Golomb ruler can be taken to represent marks 
('dots') on a physical ruler occurring at integer differences 
from each other The fact that the differences are all distinct 
implies that if a Golomb ruler is placed on top of a second, 
identical, ruler then at most one mark from the upper ruler 
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will coincide with a mark from the lower ruler, unless they 
are exactly superimposed. Golomb rulers arise in the literature 
from both theoretical and practical aspects (see [1], [4]-[6]). 
It is well known that the largest order of a ruler of length n 
is V^ + o(v^), see [4], [5], [7]. 

There are various generalizations of one-dimensional rulers 
into two-dimensional arrays. One of the most general was 
given by Robinson [8]. A two-dimensional ruler is an n x /c 
array with m dots such that all ('2 ) lines connecting two dots 
in the array are distinct as vectors, i.e., any two have either 
different length or slope. These arrays were also considered 
in [9], [10]. The case where n — k was first considered sug- 
gested by Costas and investigated by Golomb and Taylor [6]. 
Costas considered the case when n = k and each row and each 
column in the array has exactly one dot [6]. These arrays have 
application to a sonar problem [6], [11] and also to radar, 
synchronization, and alignment; they are known as Costas 
arrays. Sonar sequences are another class of arrays mentioned 
in [6], where m = k and each column has exactly one dot; 
see [12]-[14]. 

Other two-dimensional generalizations of a Golomb ruler 
have been considered in the literature, but do not have direct 
connection to our current work. For the sake of completeness 
we will mention them. A radar array is an rt x fc array with 
exactly one dot per column such that there are no two lines 
connecting two disjoint pairs of dots, occupying the same 
rows, which have the same length and slope. Radar arrays 
were defined in [6] and considered in [8], [15]-[17]. Arrays 
in which all lines have distinct slopes were considered in [18]- 
[20]. Arrays in which the Euclidean distances of any pair of 
lines are distinct were considered in [21]. 

The examples above are concerned with dots in an (infinite) 
square grid that are restricted to lie in a given line segment, 
square or rectangle. More generally, we can define a set of 
dots in a grid to be a distinct difference configuration (DDC) 
if the lines connecting pairs of dots are different either in 
length or in slope. Having surveyed the known structures of 
two-dimensional patterns with distinct differences, it seems 
that the following natural question has not been investigated: 
what is the maximum number of dots that can be placed on 
a two-dimensional square grid such that all lines connecting 
two dots are different either in their length or their slope 
and the distance between any two dots is at most r? In 
other words, rather than considering the traditional rectangular 
regions of the square grid, we consider dots which lie in 
maximal anticodes of diameter r. (An anticode of diameter 
r is a set of positions in the grid such that any pair of 
positions are at distance at most r. See Section |lll] for details.) 
We will consider two notions of distance in the square grid: 



2 



Manhattan and Euclidean. We also consider distinct difference 
configurations in the hexagonal grid, using hexagonal distance 
(also known as Manhattan distance) and Euclidean distance. 

Our motivation for considering these configurations comes 
from a new application to key predistribution for wireless 
sensor networks. We considered in [22] a key predistribution 
scheme based on DDCs in general, and Costas arrays in 
particular. A DDC A with m dots was shifted over the 
two-dimensional square grid. For each shift we assigned the 
same key to the m entries of the two-dimensional grid which 
coincide with the m dots of A. In [22] we noted that a 
Costas array is a DDC, and gave examples of DDCs with 
small numbers of dots. However, the questions of finding more 
general constructions, and providing bounds on the number of 
dots in such configurations, were left open; it is these issues 
that are addressed by the results of this paper. Other properties 
of DDCs motivated by this application are considered in [23]. 

The rest of this paper is organized as follows. In Section HIl 
we describe the models on which we will consider our two- 
dimensional patterns with distinct differences. We consider 
two two-dimensional grids: the square grid and the hexagonal 
grid. In the square grid we consider the Manhattan distance 
and the Euclidean distance, while in the hexagonal grid we 
consider the hexagonal distance and the Euclidean distance. 
We define the classes of DDCs we will study, and list 
optimal examples for small parameter sizes. In Section Hill we 
explain the relation between DDCs and maximal anticodes. 
We classify the maximal anticodes when we use Manhattan 
distance and hexagonal distance. We also briefly review some 
properties of anticodes in using Euclidean distance: these 
properties will allow us to bound the size of an anticode in 
either grid when we use Euclidean distance. In Section HVl we 
present upper bounds on the number of dots in a DDC when 
we restrict the dots to lie in some simple regions ('shapes') 
in the grid. The most important shapes we consider are the 
anticodes, in particular the Lee sphere and the hexagonal 
sphere. As a consequence of our upper bound, we settle an 
old question of Golomb and Taylor [24] (on the existence of 
honeycomb arrays of arbitrarily large size) in the negative. In 
Sections |V] and |Vl] we turn our attention to constructions and 
lower bounds for the number of dots in the classes of DDC 
defined in Section [III We generalize a folding technique that 
was used by Robinson [9] to construct Golomb rectangles, 
and provide more good examples by constructing periodic 
infinite arrays that are locally DDCs. Our constructions are 
asymptotically optimal in the case of the square grid and 
Manhattan distance. 

II. Grids, Distances, and DDCs 

We first define some new classes of two-dimensional distinct 
difference configurations. We believe that the definitions are 
very natural and are of theoretical interest, independently of 
the application we have in mind. We will consider the square 
grid and the hexagonal grid as our surface. We start with a 
short definition of the two models. Before the formal definition 
we emphasize that we define a point to be the point in 
column i and row j of either a coordinate system or a DDC. 
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Fig. 1. The hexagonal model translation 

Hence, rows are indexed from bottom to top in increasing 
order; columns are indexed from left to right in increasing 
order (So this is the usual convention for a Cartesian co- 
ordinate system, but is not the standard way of indexing the 
entries of a finite array.) 

A. The two models 

The first model is called the square model. In this model, 
a point («, j) G 1? has the following four neighbors when we 
consider the model as a connected graph: 

We can think of the points in 1? as being the centres 
of a tiling of the plane by unit squares, with two centres 
being adjacent exactly when their squares share an edge. The 
distance c?((«i,.7i), («2, ^2)) between two points and 
(*2, J2) in this model is the Manhattan distance defined by 

d{{ii,3i), (j2,j2)) = |«2 - ii\ + \h - jil- 

The second model is called the hexagonal model. Instead 
of the square grid, we define the following graph. We start by 
tiling the plane with regular hexagons whose sides have 
length l/\/3 (so that the centres of hexagons that share an 
edge are at distance 1). The vertices of the graph are the centre 
points of the hexagons. We connect two vertices if and only if 
their respective hexagons share an edge. This way, each vertex 
has exactly six neighbouring vertices. 

We will often use an isomorphic representation of the 
hexagonal model which will be of importance in the sequel. 
This representation has 1? as the set of vertices. Each point 
[x^y] e 1? has the following neighboring vertices, 

{(a: + a,y + fe) I a,6e {-1, 0, 1}, a + 7^ 0}. 

It may be shown that the two models are isomorphic by using 
the mapping ^ : ^ M^, which is defined by ^(x, y) = 
(a; + ^). The effect of the mapping on the neighbor set is 
shown in Fig. [T] From now on, slightly changing notation, we 
will also refer to this representation as the hexagonal model. 
Using this new notation the neighbors of point are 

The hexagonal distance d{x, y) between two points x and 
y in the hexagonal grid is the smallest r such that there exists 
a path with r + 1 points x ~ pi,p2, . . . ,Pr+i — y, where Pi 
and pi+i are adjacent points in the hexagonal grid. 
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B. Distinct difference configurations 

We will now define our basic notation for the DDCs we 
will focus on. 

Deflnition 1 . A Euclidean square distinct difference configu- 
ration DD(m, r) is a set ofm dots placed in a square grid such 
that the following two properties are satisfied: 

1) Any two of the dots in the configuration are at Euclidean 
distance at most r apart. 

2) All the (™) differences between pairs of dots are distinct 
either in length or in slope. 

We will also study three more classes of DDCs: A 
square distinct difference configuration DD(TO,r) is defined 
by replacing 'Euclidean distance' by 'Manhattan distance' 
in Definition [Tl a Euclidean hexagonal distinct difference 
configuration DD*(m, r) is defined by replacing 'square grid' 
by 'hexagonal grid' in Definition [T] a hexagonal distinct 
difference configuration DD (m, r) is defined by replacing 
'square grid' by 'hexagonal grid', and 'Euclidean distance' 
by 'hexagonal distance' in Definition [T] 

In the application in [22], dots in the DDC are associated 
with sensor nodes, and their position in the square or hexag- 
onal grid corresponds to a sensor's position. The parameter r 
corresponds to a sensor's wireless communication rage. So 
the most relevant distance measure for the application we 
have in mind is the Euclidean distance. Moreover, as the best 
packing of circles on a surface is by arranging the circles in a 
hexagonal grid (see [25]), the hexagonal model may be often 
be better from a practical point of view. But the Manhattan 
and hexagonal distances are combinatorially natural measures 
to consider, and our results for these distance measures are 
sharper. Note that Manhattan and hexagonal distance are both 
reasonable approximations to Euclidean distance (hexagonal 
distance being the better approximation). Indeed, since the 
distinct differences property does not depend on the distance 
measure used, it is not difficult to show that a DD(to, r) 
is a DD(m, 7'), and a DD(?Ti,r) is a DD(m, [\/2r]). Sim- 
ilarly a DD (m, r) is a DD* (m,r), and a DD* (m,r) is an 
DD*(m, r(2/V3)r]). 

C. Small parameters 

For small values of r, we used a backtrack search to 
exhaustively find a DD(7Ti,r) with m as large as possible. 
The search shows that for r = 2, 3, the largest such m are 3 
and 4 respectively, and for 4 ^ r ^ 11 the largest possible ra 
is r + 2. Fig. |2] contains examples of configurations meeting 
those bounds. 

Similarly, we found the best configurations DD (?Ti, r) in 
the hexagonal grid (see Figure [3]i for 2 ^ r ^ 10. 

III. Anticodes and DDCs 

In this section we will show a trivial connection between 
DDCs and maximal anticodes. This leads to a short investi- 
gation of maximal anticodes in the square and the hexagonal 
grids. We find all maximal anticodes in these two models under 
the Manhattan and hexagonal distance measures respectively. 




Fig. 2. Square distinct difference configurations with the largest number of 
dots possible for r = 2, 3, . . . , 11. 




Fig. 3. Hexagonal distinct difference configurations with the largest number 
of dots possible for r = 2, 3, . . . , 10. 



An anticode of diameter r in the two-dimensional grid 
(square or hexagonal) is a set S of points such that for 
each pair of points x,y E S we have d{x, y) ^ r, where 
the distance can be Manhattan, hexagonal, or Euclidean. An 
anticode S of diameter r is said to be optimal if there is 
no anticode S' of diameter r such that > |iS|. An 

anticode S of diameter r is said to be maximal if {x} U S 
has diameter greater than r for any x ^ S. Anticodes are 
important structures in various aspects of coding theory and 
extremal combinatorics [26]-[32]. 

The following two results provide an obvious connection 
between DDCs and anticodes. 

Lemma 1. Any anticode S of diameter r is contained in a 
maximal anticode S' of diameter r. 

Corollary!. A DD (m,r) is contained in a maximal anticode 
of (Euclidean) diameter r. The same statement holds for a 
DD(m, r), DD* (m,r) or DD {■m,r) when the appropriate 
distance measure is used. 
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(a) (b) (c) 

Fig. 4. Maximal anticodes in the square grid 

A. Maximal anticodes in the square grid 

We start by defining three shapes in the square grid. We 
will prove that these shapes are the only maximal anticodes 
in the square grid when we use Manhattan distance. 

A Lee sphere of radius R is the shape in the square model 
which consists of one point as centre and all positions of 
Manhattan distance at most R from this centre. The area of 
this Lee sphere is 2E? + 2_R+ 1. For the seminal paper on Lee 
spheres see [33]. Figure illustrates a Lee sphere of radius 
2. 

A bicentred Lee sphere of radius R is the shape in the square 
model which consists of two centre points (a 2 x 1 or an 1 x 2 
rectangle) and all positions of Manhattan distance at most R 
from at least one point of this centre. The area of this bicentred 
Lee sphere is 2R? + AR + 2. These shapes were used for two- 
dimensional burst-correction in [34]. Figure illustrates a 
bicentred Lee sphere of radius 2. 

A quadricentred Lee sphere of radius R is the shape in 
the square model which consists of four centre points (a 
2x2 square) and all positions of Manhattan distance at most 
R — 1 from at least one point of this centre. The area of 
this quadricentred Lee sphere is 2R? + 2R. These shapes 
were defined using the name 'generalized Lee sphere' in [35]. 
Figure 2}; illustrates a quadricentred Lee sphere of radius 3. 

Theorem 3. 

• For even r there are two different types of maximal anti- 
codes of diameter r in square grid: the Lee sphere of radius 
I and the quadricentred Lee sphere of radius |. 

• For odd r there is exactly one type of maximal anticode of 
diameter r in the square grid: the bicentred Lee sphere of 
radius . 

Proof: Let ^ be a maximal anticode of diameter r in the 
square grid. 

Assume first that r is even, so r = 2p. We will embed A 
in the two-dimensional square grid in such a way that there is 
position in A on the line y ~ x, but no position below it. The 
Manhattan distance between a point on the line y ^ x + 2p 
and a point on the line = x is at least 2p and hence A is 
bounded by the lines y ^ x and y = x + 2p. Similarly, without 
loss of generality we can assume that there is a position in A 
on the line y = —x or y = —x + 1 and no position below this 
line, so A is bounded by the lines y = —x and y = —x + 2p, 
or by the lines y = —x + 1 and y = —x + 2p+l. These four 
lines define a Lee sphere of radius p or a quadricentred Lee 
sphere of radius p. 

Now, assume that r is odd, so r = 2p + 1. We will embed 
A in the two-dimensional square grid in a way that a position 



of A lies on the line y = x, but no position lies below it. 
The Manhattan distance between a point on the line y = x + 
2/9+1 and a point on the line y = a; is at least 2p + 1 and 
hence A is bounded by the lines y = x and y = x + 2p + 1. 
Similarly, without loss of generality we can assume that there 
is a position of A on the line y = —x or y = —x + 1 and 
no position below this line, and so A is bounded by the lines 
y = —X and y = —x + 2p + 1 or by the lines y — ~x + 1 
and y = —x + 2p + 2. In either case, these four lines define a 
bicentred Lee sphere of radius p. ■ 
Finally, the following theorem is interesting from a theoret- 
ical point of view. 

Theorem 4. There exists a DD(m, r) for which the only max- 
imal anticode of diameter r containing it is a Lee sphere 
(bicentred Lee sphere, quadricentred Lee sphere) of diameter 
r. 

Proof: We provide the configurations that are needed. All 
the claims in the proof below are readily verified and left to 
the reader. 

When r is odd, we may take two points on the same 
horizontal line such that d{x,y) = r: this pair of points is 
in a bicentred Lee sphere of radius . When r is even, the 
same example is contained in a Lee sphere of radius r/2, but 
is not contained in a quadricentred Lee sphere of radius r/2. 

Let r be even, and set R = r/2. The points {0,R — 1), 
{0,R), {2R - 2, 0), {2R - 2, 2i? - 1), and (2i? - 1, R) form 
DD(5, 2R). This set of points is not contained in a Lee sphere 
of radius R, but is contained in a quadricentred Lee sphere of 
radius R. ■ 

B. Maximal anticodes in the hexagonal grid 

Theorems. There are exactly f^-j-'-] different types of maxi- 
mal anticodes of diameter r in the hexagonal grid, namely the 
anticodes Ao,Ai, . . . , A^r-i^ defined in the proof below. 

Proof: We consider the translation of the hexagonal grid 
into the square grid. By shifting it appropriately, any maximal 
anticode A of diameter r can be located inside an (r+l) x (r+ 
1) square B with corners at (0, 0), (0, r), (r, 0), and (r, r). Let 
i be defined by the property that the lines y ^ x — i contains 
a point of A, but no point of A lies below this line. 

We claim that i ^ 0. To see this, assume for a contradiction 
that i < 0. Then A is contained in the region of B bounded 
by the lines y = x — i, y ^ r, and a; = 0. But the point 
(0, r + 1) outside B is within distance r from all the points of 
this region, contradicting the fact that ^ is a maximal anticode. 
Thus, i ^ and our claim follows. 

All the points on the line y ~ x — i that are inside B 
are within hexagonal distance r from all points on the lines 
y = x — i+j,0^j^r, that lie inside B. All the points on 
the line y = x — i inside B have hexagonal distance greater 
than r from all the points on the line y = x — i + r+1. Hence, 
as A is maximal, A consists of all the points bounded by the 
lines y = x — i and y = x — i + r inside B. It is easy to verify 
that each one of the r + l anticodes Ai defined in this way is 
a maximal anticode. One can readily verify that Ai and Ar^i 



Fig. 5. Hexagonal sphere of radius 2 



(a) (b) (c) 

Fig. 6. Anticodes in the Euclidean distance 



are equivalent anticodes, since Ar-i is obtained by rotating 
Ai by 180 degrees. So the theorem follows. ■ 

Theorem 6. Let i be fixed, where ^ i ^ [■^^l • There exists 
a DD (to, r) for which the only maximal anticode of diameter 
r containing it is of the form Ai . 

Proof: Again, we provide the configurations, and leave 
the verification of the details to the reader 

For each i, 1 ^ i ^ Ai has six corner points. If 

we assign a dot to each corner point then we will obtain a 
DDC which cannot be inscribed in another maximal anticode. 
When r = 2R these six points do not define a DDC in Ab- 
In this case we assign seven dots to Ar as follows. In four 
consecutive corner points we assign a dot; in the next corner 
point we assign two dots in the adjacent points on the boundary 
of Ar; in the last corner point we assign a dot in the adjacent 
point on the boundary of Ar towards the first corner point. 
When i ~ Q, Ai is a triangle and has three corner points. If 
we assign a dot to each of these corner points we will obtain a 
DDC which cannot be inscribed in another maximal anticode. 

■ 

We now consider some basic properties of these \^-^] 
anticodes. First, the number of grid points in Ai is (r + 

1)2 _ _ ^ (i±i^ + _ The 

smallest anticode is an isosceles right triangle with base 
and height of length r + 1 containing points. The 

largest anticode is the hexagonal sphere A^r^^ of radius r/2. 

The hexagonal sphere contains ^^'^'^^ points when r is odd, 
and contains 3''^+6'-+4 ^ 3 (0^ _|_ 3 (0 + 1 points when r 
is even. The hexagonal sphere of radius R is the shape in 
the hexagonal model which consists of a centre point and all 
positions in hexagonal distance at most R from this centre 
(Fig.EJ. 

C. Maximal anticodes with Euclidean distance 

It seems much more difficult to classify the maximal an- 
ticodes in the square and hexagonal grids when we use Eu- 
clidean distance. Note that the representation of the hexagonal 
grid in the square grid does not preserve Euclidean distances, 
and so we cannot use the map ^. We expect that the overall 
shape of a maximal anticode in both models should be similar, 
since a maximal anticode in both models is just the intersection 
of a maximal anticode in with the centres of our squares or 
hexagons respectively. But the iocaF structure of an anticode 
will be different; for example, in the hexagonal grid we can 
have three dots that are pairwise at distance r, but this is not 
possible in the square grid. 



Because maximal anticodes in determine the shape 
of maximal anticodes in the square or hexagonal models, 
we conclude this section with a brief description of such 
anticodes. 

An anticode is confined to the area as depicted in Figure |6^, 
where dots are two elements in the anticode at distance r. 
The most obvious maximal anticode is a circle of diameter r 
depicted in Figure ^p. Another maximal anticode is depicted 
in Figure |6};, and is constructed by taking three dots at the 
vertices of an equilateral triangle of side r, and intersecting the 
circles of radius r about these dots. Between the 'triangular' 
anticode and the circle there are infinitely many other max- 
imal anticodes. We will need the following 'isoperimetrical' 
theorem; see Littlewood [36, Page 32] for a proof. 

Theorem 7. Let Abe a region ofM? of diameter r and area a. 
Then a < (7r/4)r2. 

We remark that the example of a circle of diameter r shows 
that the bound of this theorem is tight. 

IV. Upper bounds on the Number of Dots 

In this section we will provide asymptotic upper bounds 
on the number of dots that can be contained in a DDC, 
using a technique due to Erdos and Turan [18], [37]. We 
start by considering upper bounds on the number m of dots 
in a DD(TO,r) and a DD {m,r), and then consider upper 
bounds in a DD(m, r) and DD*(TO,r). The results for small 
parameters in Section |ll] might suggest that a DD(7Ti,r) 
can always contain r + 2 dots: our result (Theorem |9]l that 
TO ^ -^r + o{r) surprised us. Our techniques easily generalise 
to DDCs where we restrict the dots to lie in various shapes in 
the grid not necessarily related to distance measures: we end 
the section with a brief discussion of this general situation. 

A. Manhattan and hexagonal distances 

Lemma 8. Let r bea non-negative integer Let A be an anticode 
of Manhattan diameter r in the square grid. Let £ be a positive 
integer such that I ^ r, and let w be the number of Lee spheres 
of radius £ that intersect A non-trivially. Then w ^ ^{r+2£)'^ + 
0{r). 

Proof: Let A' be the set of centres of the Lee spheres 
we are considering, so w = \A'\. We claim that A' is an 
anticode of diameter at most r + 2£. To see this, let c, c' e A'. 
Since the spheres of radius £ about c and c' intersect A non- 
trivially, there exist elements a,a' G A such that d{c, a) ^ £ 
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and d{c', a') ^ £. But then 

d(c, c') < d{c, a) + d{a, a') + d{a' , c') £ + r + £ = r + 2£, 

and so our claim follows. 

Let A" be a maximal anticode of diameter r + 2£ containing 
A'. Theorem |3] implies that A" is a Lee sphere, bicentred 
Lee sphere or quadricentred Lee sphere of radius R, where 
R = [{r + 2£)/2l. In all three cases, \A"\ = 2R^ + 0{R) = 
\[r + 2£f + 0{r). Since 



l^'l 



the lemma follows. 



Theorem 9. If a DD(m, r) exists, then 

^^ + (3/24/3)^2/3^0(^1/3)^ 

Proof: We begin by giving a simple argument that leads 
to a linear bound on 77i in terms of r, with an inferior leading 
term to the bound in the statement of the theorem. There are 
2r^ + 2r non-zero vectors of Manhattan length r or less, where 
a vector is a line with direction which connects two points. 
The distinct difference property implies that each such vector 
arises at most once as the vector difference of a pair of dots. 
Since a configuration of m dots gives rise to m{m— 1) vector 
differences, we find that 

m.{m - 1) ^ 2r^ + 2r. 

In particular, we see that m ^ \/2r + o(r) = 0{r). 

We now estabUsh the bound of the theorem. Since all the 
dots are at distance at most r, we see that all dots are contained 
in a fixed anticode A of the square grid of diameter r. Set 
£ = [ar^/'^J, where we will choose the constant a later so 
as to optimize our bound. We cover A with all the 'small' 
Lee spheres of radius £ that intersect A nontrivially. Every 
point of A is contained in exactly a small Lee spheres, where 
a — 2£^ + 2^+1. Moreover, by Lemma [S) we have used w 
small Lee spheres, where w ^ ^(^ + 2£)^ + 0(r). 

Let mi be the number of dots in the ith small Lee sphere. 
Let /i be the mean of the integers m^. Since every dot is 
contained in exactly a small Lee spheres, /x = am/w. We 
aim to show that 



^(/i^ — /i) ^ mi{mi — 1) ^ a{a — 1). 



(1) 



The first inequality in ([T]i follows from expanding the non- 
negative sum so it remains to show the second 
inequality. 

The sum J27=i ^ii^i ~ 1) counts the number of pairs 
(£, d) where £ is a small Lee sphere and d is a vector 
difference between two dots in £. Every difference d arises 
from a unique ordered pair of dots in A, since the dots form 
a distinct difference configuration. Thus 



E 



■mi{m.i - 1) s$ ^ k{d), 
d 



where we sum over all non-zero vector differences d and where 
k{d) is the number of Lee spheres of radius £ that contain any 



fixed pair of dots with vector difference d. If we assume that 
the first element of the pair of dots with vector difference d 
lies at the origin, we see that 

^ k{d) = a{a - 1) 

d 

since there are exactly a Lee spheres of radius £ containing the 
origin, and each such sphere contributes 1 to k{d) for exactly 
a — 1 values of d. Thus we have established 

Now, the inequality ^ together with the fact that = 
am/w imply that 



and so 



(fi — l)m ^ a — 1 ^ a, 



m 5j w ( 1 + 



(2) 



By Lemma [HI 



and we have that 

+ (m/a)) = 1 + m/{2a) + 0{{m/af). 

Since m = 0{r) and a ^ 2a^7'4/3^ these two inequalities 
combine with (|2]i to show that 



m ^ 



(3) 



Since m = 0{r), this inequality implies that m ^ -^r + 
0(r^/'^). Combining this tighter bound with ([3]) we find that 



m ^ — 1 
V2 



2a 



4V2a2 



,,-1/3 



0(r 



-2/3- 



The expression 2a + \ / {4,^/2,0^) is minimized when a = 
2-5/6 jjjg value 3/2^/^, so choosing this value for a we 
deduce that 

m ^^r(l + 5^r-i/3 + 0(^-2/3)) 

= 72^+2^^'/'+0(^'^')' 

as required. ■ 
We now look at the hexagonal grid. 

Lemma 10. Let r be a non-negative integer. Let A be an 
anticode of hexagonal diameter r in the hexagonal grid. Let £ 
be a positive integer such that £ ^ r, and let w be the number 
of hexagonal spheres of radius £ that intersect A non-trivially. 
Thenw |(r + 2£)2 + 0(r). 

Proof: The set of centres of the hexagonal spheres of 
radius £ that have non-trivial intersection with A clearly form 
an anticode of diameter at most r + 2£. Therefore the number 
w of such spheres is bounded by the maximal size of such 
an anticode. The results on the maximal anticodes in the 
hexagonal metric in Section [III] imply that 

w < jiHr + 2£f + 6(r + 2£) + 4) 



Kr + 2£)2 + 0(r), 



as required. 
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Theorem 11. If a DD (m, r) exists, then 

m^^r+ (34/32-5/3)^2/3 ^ o{r^/'). 

Proof: The dots in a DD (m, r) form an anticode of 
diameter r. Let ^ = [2^^/'^3^^/^r^/'^J . We may cover these 
dots with the w hexagonal spheres of radius £ that contain 
one or more of these dots. By Lemma [TO] we have that 
w ^ |(r + 2^)2 + 0(r). 

Using the fact that a hexagonal sphere of radius £ contains a 
points in the hexagonal grid, where a = + 3^+1, we may 
argue exactly as in Theorem |9] to produce the bound (|2]i. There 
are O(r^) vectors in the hexagonal grid of hexagonal length 
r or less, so the argument in the first paragraph of Theorem |9] 
shows that m = 0{r). Since m/a = 0{r^'^^'^), the bound (|2]i 
implies that 

m ^ + 0(r2/3) = + 0[t'''^). 

This bound on m implies that m/a = 2^/'^3^^/^r^^/'' + 
(9(7.-2/3)^ and so applying (|2]i once more we obtain the bound 
of the theorem, as required. ■ 
One consequence of Theorem [TT| is an answer to the ninth 
question asked by Golomb and Taylor [24]: a honeycomb army 
is a DDC in the hexagonal grid whose dots, when represented 
in the square grid, form an m x m Costas array whose 
dots lie in m consecutive 'North-East' diagonals. Honeycomb 
arrays are the natural hexagonal analogue of Costas arrays. Do 
honeycomb arrays exist for infinitely many ml The conjecture 
in [24] is that the answer is YES. However, the answer is in 
fact NO, as the following corollary to Theorem [TTI shows. 

Corollary 12. Honeycomb arrays exist for only a finite number 
of values ofm. 

Proof: The dots in a honeycomb array are contained 
in an anticode of diameter at most m — 1 (using hexagonal 
distance). Hence a honeycomb array is a DD [m, m— 1). But 
TheoremfTT] shows that m s$ ^m + Oim^/^). Since ^ < 1, 
no honeycomb array exists when m is sufficiently large. ■ 
In fact, numerical computations indicate that no honeycomb 
arrays exist for in ^ 1289: for m in this range, there is a 
suitable choice of £ such that a honeycomb array violates the 
inequality Q. 

B. Euclidean distance 

We now turn our attention to Euclidean distance. Our first 
lemma is closely related to Gauss's circle problem: 

Lemma 13. Let £ be a positive integer, and let S be a (Eu- 
clidean) circle of radius £ in the plane. Then the number of 
points of the square grid contained in S is + 0{£). 

Proof: Let c be the centre of S. Let X be the set of points 
of the square grid contained in S. Define X to be the union 
of all unit squares whose centres lie in X. Clearly X has area 
\X\. The maximum distance from the centre of a unit square 
to any point in the unit square is at most l/%/2, and so X 
is contained in the circle of radius £ + {\/^/2) with centre c. 



Similarly, every point in a circle of radius £ — (l/\/2) with 
centre c is contained in X. Hence 

4£ - (l/%/2)f < \X\ sc ^{£ + (1/V2))2, 

and so the lemma follows. ■ 

Lemma 14. Let r be a non-negative integer. Let A be an 
anticode in the square grid of Euclidean diameter r. Let £ be 
a positive integer such that £ ^ r, and let w be the number of 
circles of radius £ whose centres lie in the square grid and that 
intersect A non-trivially. Then w ^ (7r/4)(r + 2£)^ + 0{r). 

Proof: As in Lemma [H it is not difficult to see that the set 
A' of centres of circles we are considering form an anticode in 
the square grid of diameter at most r + 2£. Note that w — \A'\. 
Let X be the union of the unit squares whose centres lie in A', 
so X has area w. The maximum distance between the centre 
of a unit square and any other point in this square is I/V2, and 
so X is an anticode in of diameter at most r + 2£+{l/ \/2) . 
Hence, by Theorem|71 w ^ (7r/4)(r + 2^+ (1/V2))2 and the 
lemma follows. ■ 

Theorem 15. If a DD(to, r) exists, then 

Proof: The proof is essentially the same as the proof 

of Theorem (TT] using Lemma [13] to bound the number a of 

points in a sphere of radius £, and using Lemma [14] instead 

of Lemma [TOl The bound of the theorem is obtained if we set 
£= [l/(22/37ri/6)7.2/3j, I 

Lemma 16. Let £ be a positive integer, and let S be a (Eu- 
clidean) circle of radius £ in the plane. Then the number of 
points of the hexagonal grid contained in S is (27r/y3)£2 + 
0{£). 

Proof: The proof of the lemma is essentially the same as 
the proof of Lemma [T3] The hexagons whose centres form the 
hexagonal grid have area V3/2, and the maximum distance 
of the centre of a hexagon to any point in the hexagon is 
1/V3- Define X to be the set of points of the hexagonal grid 
contained in S, and let X be the union of all hexagons in our 
grid whose centres lie in X. Clearly X has area (\/3/2)|X|. 
The argument of Lemma [13] shows that 

n{£ - (l/VS))^ ^ iV3/2)\X\ ^ n{£ + (l/VS))^, 

and so the lemma follows. ■ 

Lemma 17. Let r be a non-negative integer. Let A be an 
anticode in the square grid of Euclidean diameter r. Let £ be 
a positive integer such that £ ^ r, and let w be the number of 
circles of radius £ whose centres lie in the hexagonal grid and 
that intersect A non-trivially. Then w ^ {n/ {2^)){r + 2£)'^ + 
0{r). 

Proof: The proof of this lemma is essentially the same as 
the proof of Lemma [14] The argument there with appropriate 
modifications shows that {V3/2)w ^ (7r/4)(r+2£+(l/V3))2 
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(where the factor of \/3/2 comes from the fact that the 
hexagons associated with our grid have area \/3/2). ■ 

Theorem 18. If a DD* (m, r ) exists, then 

Proof: The proof is the same as the proof of Theorem [TSl 
using Lemmas [16] and [17] in place of Lemmas [13] and [14] and 
defining e = L3i/i22-V6^-i/6^2/3j , 



C. More general shapes 

All the theorems above consider a maximal anticode in 
some metric, and cover this region with small circles of radius 
£. We comment (for use later) that the same techniques work 
for any 'sensible' shape that is not necessarily an anticode. 
(We just need that the number of small circles that intersect 
our shape is approximately equal to the number of grid points 
contained in the shape.) So we can prove similar theorems 
for DDCs that are restricted to lie inside regular polygons, for 
example. The maximal number of dots in such a DDC is at 
most y/s + o{y/s) when the shape contains s points of the 
grid. 

V. Periodic Two-Dimensional Configurations 

The previously known constructions for DDCs restrict dots 
to lie in a line or a rectangular region (often a square region) of 
the plane. The application described in [22] instead demands 
that the dots lie in some anticode. The most straightforward 
approach to constructing DDCs for our application is to find 
a large square or rectangular subregion of our anticode, and 
use one of these known constructions to place dots in this 
subregion. This approach provides a lower bound for m that 
is linear in r, but in fact we are able to do much better 
than this by modifying known constructions (in the case 
of Robinson's folding technique below) and by making use 
of certain periodicity properties of infinite arrays related to 
rectangular constructions. We will explain how this can be 
done in the next section. In this section we will survey some 
of the known constructions for rectangular DDCs, extend 
these constructions to infinite periodic arrays, and prove the 
properties we need for Section [VT] 

Let A he a (generally infinite) array of dots in the square 
grid, and let rj and k be positive integers. We say that A is 
doubly periodic with period (77, k) if A{i, j) = A{i + i], j) and 
^(i, j) = -4(i,j + k) for all integers i and j. We define the 
density of A to be d/{r]K), where d is the number of dots in 
any k x ?/ sub-array of A. Note that the period (jy, k) will not 
be unique, but that the density of A does not depend on the 
period we choose. We say that a doubly periodic array A of 
dots is a doubly periodic nx k DDC if every nx k sub-array 
of is a DDC. See [14], [38], [39] for some information 
on doubly periodic arrays in this context. We aim to present 
several constructions of doubly periodic DDCs of high density. 



A. Constructions from Costas Arrays 

A Costas array of order n is an ?i x 71 permutation array 
which is also a DDC. Essentially two constructions for Costas 
arrays are known, and both give rise to doubly periodic DDCs. 

The Periodic Welch Construction: 

Let a be a primitive root modulo a prime p and let A be the 
square grid. For any integers i and j, there is a dot in Aii,j) 
if and only if a* = j mod p. 

The following theorem is easy to prove. A proof which also 
mentions some other properties of the construction is given 
in [23]. 

Theorem 19. Let A be the array of dots from the Periodic 
Welch Construction. Then A is a doubly periodic p x {p — 1) 
DDC with period [p — 1, p) and density 1/p. 

Indeed, it is not difficult to show that each p x [p ~ 1) sub- 
array is a DDC with p — 1 dots: a dot in each column and 
exactly one empty row. The {p — 1) x {p — 1) sub-array with 
lower left corner at .4(1, 1) is a Costas array. 

The Periodic Golomb Construction: 

Let a and /3 be two primitive elements in GF{q), where q 
is a prime power. For any integers i and j, there is a dot in 
AiiJ) if and only if a' + f3^ = 1. 

The following theorem is proved similarly to the proof 
in [23], [40]. 

Theorem 20. Let A be the array of dots from the Periodic 
Golomb Construction. Then A is a doubly periodic (q— 1) x (q— 
1) DDC with period {q-l,q-l) and density {q-2)/{q- 1)^. 

Indeed, each (q — 1) x (q — 1) sub-array of ^ is a DDC with 
q — 2 dots; exactly one row and one column are empty. The 
(g — 2) X (g — 2) sub-array with lower left corner at ^(1, 1) 
is a Costas array. 

If we take a = /3 in the Golomb construction, then the 
construction is known as the Lempel Construction. There are 
various variants for these two constructions resulting in Costas 
arrays with orders slightly smaller (by 1, 2, 3, or 4) or larger 
by one than the orders of these two constructions (see [24], 
[41]). These are of less interest in our discussion, as they do 
not extend to doubly periodic arrays in an obvious way. 

B. Constructions from Golomb rectangles 

A Golomb rectangle is an n x fc DDC with m dots; Costas 
arrays are a special case. Apart from constructions of spe- 
cial cases, there is essentially one other general construction 
known, the folded rulers construction due to Robinson [9]. 
Folded Ruler Construction: 

Let S = {ai, a2, • • • , a™} C {0, 1, ... , n} he a Golomb 
ruler of length n. Let £ and k he integers such that £-k ^ n+1. 
Define A to be the £ x k array where A{i,j), ^ i ^ k — 1, 
^ J < ^ — 1, has a dot if and only if i- £ + j = at for some 
t. 

Theorem 21. The array A of the Folded Ruler Construction is 
antxk Golomb rectangle. 
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We now show how to adapt the Folded Ruler Construction 
to obtain a doubly periodic £x k DDC. We require a stronger 
object than a Golomb ruler as the basis for our folding 
construction, defined as follows. 

Definition 2. Let A be an abelian group, and let D = 
{ai, a2, ■ ■ ■ , am} C A be a sequence of m distinct elements 
of A. We say that D is a _B2-sequence over A if all the sums 
a.i-^ + with 1 ^ ii ^ 12 ^ TO are distinct. 

For a survey on _B2-sequences and their generalizations the 
reader is referred to [42]. The following lemma is well known 
and can be readily verified. 

Lemma 22. A subset D = {ai, 02, ... , a,„} C A is a B2- 
sequence over A if and only if all the differences — with 
l^ii7^Z2^m are distinct in A. 

So, in particular, a Golomb ruler is exactly a i32 -sequence 
over Z. Note that a i?2-sequence {ai, a2, . . . , over Z„ 
produces a Golomb ruler {fei, 62, ■ • • , whenever the 6; are 
integers such that a; = bi mod n. Also note that if D is a _B2- 
sequence over Z„ and a € Z„, then so is the shift a + D ~ 
{a + d : d E D}. The following theorem, due to Bose [43], 
shows that large i32-sequences over Z„ exist for many values 
of n. 

Theorem 23. Let q be a prime power Then there exists a B2- 
sequence ai,a2, . . . , am over Z„ where n ^ — 1 and m = q. 

The Doubly Periodic Folding Construction: 

Let n be a positive integer and D = {ai, 02, ... , a,„} be a 
i32-sequence in Z„. Let £ and k be integers such that £-k ^ n. 
Let A be the square grid. For any integers i and j, there is a 
dot in A{i,j) if and only if at = i ■ £ + j mod n for some t. 

Theorem 24. Let A be the array of the Doubly Periodic Folding 
Construction. Then A is a doubly periodic £x k DDC of period 

i gcd(n,e) ' ^) d density m/n. 

Proof: Let /(x, y) ~ x ■ £ + y. The period of A follows 
from the observation that for each two integers a and /3 we 
have f{i,j) = f{i + a-^^^(^J + Pn) = i-£ + j mod n. The 
density of A is m/n follows since there are exactly to dots 
in any n consecutive positions in any column. 

Let S he an £ X k sub-array, whose lower left-hand corner 
is at A{i,j). An alternative construction of the dots in S is as 
follows. Take the shift {i-£ + j) + D of D, which is also a B2- 
sequence in Z„. Let D' be the corresponding Golomb ruler 
in {0, 1, . . . , 71 — 1}, so a e if and only if a = 6 mod n, 
where b £ {i ■ £ + j) + D. Then form dots in S by using the 
Folded Ruler Construction. Hence, by Theorem [21] the dots 
in S form a DDC and so the theorem follows. ■ 

The following slightly different construction also produces 
doubly periodic Golomb rectangles. 

The Chinese Remainder Theorem Construction: 

Let n be a positive integer and let D = {oi, 02, . . . , a,„} 
be a i?2 -sequence in Z„. Let n = £ • fc be any factorization of 



n such that gcd(^, k) = 1. For any two integers i and j we 
place a dot in A{i,j), if and only if at ~ {i ■ £ + j ■ k) mod n 
for some t. 

Theorem 25. Let A be the array constructed by the Chinese 
Remainder Theorem construction. Then A is a doubly periodic 
£ X k DDC of period {k, £) and density m/n. Moreover, every 
£ X k sub-array of A contains exactly m dots. 

Proof: Let f{x, y) = x ■ £ + y ■ k. For any two integers a 
and (3 we have f{i,j) = f{i + ak,j + p£) = i-£+j-k mod n. 
So the definition of A impUes that A is doubly periodic with 
period {k,£). 

Since £ and k are relatively primes, it follows (from the 
Chinese Remainder Theorem) that each integer s in the range 
^ s ^ £-k—l, has a unique representation as s = d-£ + e-k, 
where O^d^fc — 1, 0^e^£ — 1. Hence every £ x k sub- 
array of A has m dots corresponding to the m elements of the 
i32-sequence D. In particular, this implies that A has density 
m/n. 

Assume for a contradiction that there exists an £ x fc sub- 
array 5 of that is not a DDC. Suppose that the lower left- 
hand corner of S is at A{i,j). The distribution of dots in 
S is the same as the distribution in the sub-array with lower 
left-hand corner the origin once we replace D by the shift 
{i£ + j£) + D. So, without loss of generality, we may assume 
that the lower left-hand corner of iS lies at the origin. As the 
distinct difference property fails to be satisfied, there are four 
positions with dots in A of the form: 

A{ii,ji) A{ii+dji+e) 
-4(i2,j2) A{i2+d,j2+e) 

where + d, 12, i2 + d E {0, 1, . . . , fc — 1} and + 
e, j2, j2 + e e {0,1, ...,£— 1}. By the definition of A we 
have 

iii+d)£+{ji + e)k - {ii£ + jik) = d-£ + e-k 
{i2 + d)£ + (^2 + e)k - {i2£ + 32k) = d-£ + e-k 

Since by Lemma |22] each nonzero residue s modulo n has at 
most one representation as a difference from two elements of 
D, it follows that the pairs 

{A{ii,ji),A{ii+d,ji+e)} 

{Ail2,j2),Ail2+d,j2+e)} 

are identical, and the theorem follows. ■ 

VI. Lower bounds 

A. Manhattan distance 

In this section we will prove that there exists a DD(TO,r) 
with ^ — o(r) dots: this attains asymptotically the upper 
bound of Theorem |9] We will see that this construction is 
actually using folding in a slightly different way. We further 
show that we can construct a doubly periodic array in which 
each Lee sphere of diameter r is a DDC with + o(r) dots. 

The LeeDD Construction: 

Let r be an integer, and define R = [|J. Let D = 
{ai, 02, . . . , a^} be a ruler of length n. Define f{i,j) = 
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iR + j{R+l) + R'^ + R. Let Ahe the Lee sphere of radius R 
centred at (0,0), so A has the entry A{i,j) if \i\ + \j\ R. 
We place a dot in A{i,j) if and only if f{i,j) € D. 

Theorem 26. The Lee sphere A of the LeeDD Construction is 

aDD(7n,r), wherem^ \D n {0,1, . . . ,2R'^ + 2R}\. 

Proof: We first note that if |i| + |j| ^ i? then the smallest 
value that the function / takes is and the largest value is 
2R^ + 2R. Next, we claim that if and (12,^2) are two 



distinct points such that + \ji 



IJ2I ^ R then 



fiiiji) 7^ /(«2,i2)- Assume the contrary, that /(ii, ji) = 
/(«2,j2). SoiiR + n{R+l) + R^ + R = i2R + j2iR+l) + 
+ and therefore (12- = (ji-.?2)(^+l)- If «i = «2, 
then ji = j2 which contradicts our assumption that and 
{12,32) are distinct. So we may assume that ii 7^ 12- Similarly, 
we may assume that ji ^ j2- The equality (12 — ii)R — 
(ji ^ J2)(-R + 1) now implies that R+1 divides \i2 — ii\ and 
i? divides |j2^Ji|- This implies that 1*2 + b2 ^ji| > 2i?, 
but 



1*2 - ill + 



IJ2-J1I < Nl| + |jl| + N2| + b2K2i?, 



and so we have a contradiction. Thus, f{ii,ji) 7^ /(i2,j2)- 
This implies that each one of the integers between and 
2R^ + 2R is the image of exactly one pair In partic- 

ular, the number m of dots in the configuration is exactly 
|Dn{0,l,...,2i?2 + 2i?}|. 

Since ^ is a Lee sphere of radius R, it follows that the 
Manhattan distance between any two points is at most 2R ^ r. 
Now, assume for a contradiction that A is not a DD(7ti, r), so 
there exist four positions with dots in A as follows: 



■A{i2,32) 



A{ii 
A{i2 



d, ji 
d, j2 



e) 
e) 



By definition we have that 

f{hJi)J{ii + d,ji +e),f [12,32), f{i2 + d,j2 + e) e D. 

But then /{h + d,ji + e) - /(h, ji) = /(h + d,ji + e) - 
fihiji) = dR + e{R + 1), contradicting the fact that D is a 
ruler 

Thus, the Lee sphere A of the LeeDD Construction is a 
DD(m, r). ■ 

Corollary 27 There exists a DD (m,r) in which to = — 
o(r). 

Proof: Define R = [r/2j and let n = 2R'^+2R+1. There 
exists a ruler of length at most n containing m dots, where 
m> see [4], [5], [1]. Let DC {0,1, . . . ,n-l} 

he such a ruler The corollary now follows, by Theorem |26] 

■ 

It is worth mentioning that the LeeDD Construction is 
actually a folding of the ruler by the diagonals of the Lee 
sphere. Figure |7] illustrates why this is the case, by labelling 
the positions in a Lee sphere of radius 3 by the values of 
f{i,j) at these positions. So if we use a i32-sequence over 
Z„ instead of a ruler in the LeeDD Construction we obtain a 
doubly periodic array with nice properties: 
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Fig. 7. Folding along diagonals 



Fig. 8. A (3,5)-diagonally extended Lee sphere 



The Doubly Periodic LeeDD Construction: 

Let r he an integer, R ~ [|J, and let D = {oi, a2, . . . , a^i} 
be a i32-sequence over Z„, where n ^ 2R^ + 2R + 1. Let 
f{i,j) = iR + j{R+l) mod n. Let A he the square grid. For 
each two integers i and j, there is a dot in A{i,j) if and only 

if fit, 3) eD. 

Similarly to Theorem |26] we can prove the following result. 

Theorem 28. The array A constructed in the LeeDD Construc- 
tion is doubly periodic with period [n, n) and density n/n. The 
dots contained in any Lee sphere of radius R form a DDC 

Proof: The first statement of the theorem is obvious. The 
second statement follows as in the proof of Theorem |26] once 
we observe that / is an injection when restricted to any Lee 
sphere of radius R. ■ 
In Subsections IVI-DI and IVI-EI we will make use of an 
extension of this construction. For positive integers R and t, 
an {R,t)-diagonally extended Lee sphere is a set of positions 
in the square grid defined as follows. Let (io,.?o) € 1?, and 
define C = {(iq + k, + k) : k ^ t - 1}. Then an 
{R, t)-diagonally extended Lee sphere is the union of the Lee 
spheres of radius R with centres lying in C. (See Fig.[8]for an 
example.) An [R, t)-diagonally extended Lee sphere contains 
exactly 2R^ + t{2R+ 1) positions; the Lee sphere of radius R 
is the special case when t ~ 1. We observe that by choosing 
?7 ^ 2R^ + t{2R + 1), we can generalize the doubly periodic 
LeeDD construction by continuing folding along the diagonals 
of the rectangle. This yields the following corollary, which 
will prove useful in the construction of configurations for the 
hexagonal grid. 

Corollary 29. Let a be positive, and let n be an integer such 
that n ^ (2 + 2a)R? + aR. Consider the array A constructed 
using the doubly periodic LeeDD Construction. Then A is a 
doubly periodic array with density /i/?7. The dots contained in 
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any (R, [aR\) -diagonally extended Lee sphere form a DDC. 
There exists a family of B2 sequences so that A has density at 
least 1 / ^(2 + 2a)i?2 + o(i?2). 

Proof: To establish the final statement of the corollary, 
we choose a family of B2 sequences as follows. Let p be 
the smallest prime such that p'^ — 1 ^ (2 + 2a)R^ + ai?, 
and define n = p^ — \. By Ingham's classical result [44] on 
the gaps between primes, we have that n ^ (2 + 2a)R? + 
0[R^^/») ^ (2 + 2a)R^ + o{R^). By TheoremES there exists 
a B2 sequence over Z„ with n ~ p. Hence the density of A 



n/n = p/(p2 - 1) ^ l/v/(2 + 2a)i?2 + o(i?2)^ 




as required. 



Fig. 9. Square intersecting a circle 



B. A General Technique 

Let 5 be a shape (a set of positions) in the square grid. We 
are interested in finding large DDCs contained in S, where (for 
example) S is an anticode. This subsection presents a general 
technique for showing the existence of such DDCs, using the 
doubly periodic constructions from Section W\ 

We write (i, j) + S for the shifted copy {{i + i' ,j + f) : 
G S} of S. Let ^ be a doubly periodic array. We 
say that ^ is a doubly periodic S-DDC if the dots contained 
in every shift (i, j) +5 of S form a DDC. So the doubly 
periodic arrays constructed in Section|V]are all doubly periodic 
5-DDCs where 5 is a square or a rectangle; the arrays in 
Theorem |28] and Corollary |29] are doubly periodic 5-DDCs 
with S a Lee sphere and diagonally extended Lee sphere 
respectively. The following lemma follows in a straightforward 
way from our definitions; 

Lemma 30. Let Abe a doubly periodic S-DDC, and let S' C 
S. Then A is a doubly periodic S' -DDC. 

We will use doubly periodic DDCs to prove the existence 
of the configurations we are most interested in, using the 
following theorem. 

Theorem 31. Let S be a shape, and let Abe a doubly periodic 
S-DDC of density 5. Then there exists a set of at least \5\SW 
dots contained in S that form a DDC. 

Proof: Let the period of A be (77, k). Write mi,j for the 
number of dots of A contained in the shift j) + S of S. 
Now A is periodic, so the definition of the density of A shows 
that 



i=i i=i 

Hence the average size of the integer rriij is 5\S\, so there 
exists an integer to^'j' such that mi' ^ji ^ r<5|5|]. The TOi'j' 
dots in +S form a DDC, by our assumption on A, and 

so the appropriate shift of these dots provides a DDC in S 
with at least dots, as required. ■ 



C. Euclidean distance in the square model 

This subsection illustrates our general technique in the 
square grid using Euclidean distance. So we wish to construct 
a DD(m, r) with m as large as possible. 

Let R= [r/2], and let S be the set of points in the square 
grid that are contained in the Euclidean circle of radius R 
about the origin. We construct a DDC contained in S with 
many dots: any such configuration is clearly a DD(7Ti,r) for 
some value of m. The most straightforward approach is to find 
a large square contained in S (which will have sides of length 
approximately V2R), and then add dots within this square 
using a Costas array. This will produce a DD{m,r) where 

m = V2R - o{R) = ^r- o{r) « 0.707r. 

To motivate our better construction, we proceed as follows. We 
find a square of side n where n > \/2R that partially overlaps 
our circle: see Figure |9] The constructions of Section |V] show 
that there exist doubly periodic rix n DDCs that have density 
approximately 1/n. So Theorem [3T| shows that for any shape 
S' within the square, there exist DDCs in S' that have at least 
|5'|/?7 dots. Let S' be the intersection of our square with S. 
Defining 9 as in the diagram, some basic geometry shows that 
the area of S' is 

(7r/2) - 26* + sin 26*, 



\S'\ 



\S\ = 2R\{'k/2) -29 + sin 261). 



2 cos2 9 

Since n — 2R cos 9, Theorem [3T| shows that the density of 
dots within S' can be about 1/n = l/(2i?cos6') when n is 
large. So we can hope for at least i^iR dots, where 1^1 is the 
maximum value of 

((7r/2) - 261 + sin 261)/ cos 

on the interval ^ 6' ^ 7r/4. In fact /i w 1.61589, achieved 
when 9 w 0.41586 (and so when ?? ~ rcos9 = cr, where 
c« 0.914769). 

Theorem 32. Let 1^1 be defined as above. There exists a 
DD(m, r) in which m = ifi/2)r ~ o{r) w 0.80795r. 

Note that Theorem [15] gives an upper bound on m of the 
form TO < iV^/2)r + o(r) w 0.88623r. 
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Fig. 10. A diagonal rectangle intersecting the image of a hexagonal sphere 



Proof: Define c w 0.91477 as above. Let q be the 
smallest prime power such that q > cr. We have that 
cr < q < cr + {crf^^, by a classical result of Ingham [44] 
on the gaps between primes; so in particular q ^ cr. By 
Theorem |20l there exists a doubly periodic (q — 1) x [q — 1) 
DDC A of density {q - 2)/{q - 1)^. Let 5' be the inter- 
section between S and a Euclidean circle of radius [r/2j 
about the origin. Then ^ is a doubly periodic 5'-DDC. By 
Theorem UT\ there exists a DDC in S' with at least m dots, 
where \S'\{q—2)/{q— 1)^. But the geometric argument above 
shows that \S' \{q — 2) / {q — 1)'^ ~ (/i/2)r, and so the theorem 
follows. ■ 



D. Hexagonal distance 

By representing the hexagonal anticodes in the square grid, 
we may use TheoremlSTIto show the existence of a DD (m, ?■) 
where m is large. The method of producing lower bounds is 
essentially the same as above, but the geometrical problem 
we are solving is different, with the images under ^ of the 
maximal anticodes Ai replacing the circle, and the DDC 
contained a diagonally extended Lee sphere of Corollary |29] 
replacing the Costas array contained in a square. Here we 
consider the case of configurations contained in the hexagonal 
sphere Af(^r-i)/2y^ the cases of the other anticodes may be 
handled in a similar fashion. The problem we are solving 
is pictured in Fig. \TQ\ The figure shows the image under ^ 
of the hexagonal sphere of radius R = [r/2\ in bold; the 
square of side 2R + 1 containing this image is also shown. 
The hexagonal sphere contains a Lee sphere of radius R with 
the same centre; the region S we consider is the [R, [ai?J )- 
diagonally extended Lee sphere whose mid-point is at the 
centre of the hexagonal sphere; see Fig. [TO] Let S' be the 
intersection of S with the image of the hexagonal sphere. We 
have that \S'\ ^ R^{2 + 2a- a^) + o{R^). By Corollary |29l 
there is a doubly periodic 5-DDC of density at least 
where n — 2R^{1 + a) + o(i?^). Thus Theorem UT\ shows that 
there is a DDC contained in S' containing fiR. — o{R) dots, 
where fi is the maximum of 

2 + 2a-a^ 




Fig. 11. A (t, [aij )-diagonally extended Lee sphere is transformed into a 
rotated square (when at = (\/3- l)i + 1) 



It can be seen that /i 



1+2V7 



1.58887, achieved 



when a = — i^-^. Since S' is contained in a hexagonal sphere 
of radius R, all pairs of dots in our DDC are at hexagonal 
distance at most r. Thus we have the following theorem: 



Theorem 33. Let ^ be defined as above. There exists a 
DD*(m, r) in whicli m ifJ-/2)r - o{r) « 0.79444r. 

E. Euclidean distance in the hexagonal model 

In this subsection we will obtain a construction for a 
DD* (m,r) contained within a circle of radius R = \j/2\, 
again based on the doubly periodic LeeDD construction. We 
first observe that a diagonally extended Lee sphere in the 
square grid is transformed by into a (rotated) rectangle 
in the hexagonal grid. In particular, a (i, [(\/3 — l)t + IJ)- 
diagonally extended Lee sphere is transformed by into 
a set of hexagons whose centres all lie within a (rotated) 
square S of side \/3t (see Fig. [TTT i. Corollary |29] shows that 
there is a doubly periodic 5-DDC with density 1 / -/n, where 
n = 2Vit^ + o{t^). 

Consider (see Fig. [T2l i a circle of radius R and a square S 
of side s where s = 2i?cos 9. Since a hexagon has area a/3/2, 
the square S contains (8/-\/3)-R^ cos^ 9 + 0{R) hexagons. Let 
S' be the intersection of S with the circle of radius R. The 
calculations in Subsection I VI-CI show that 

(7r/2 - 2e' + sin26') 



\S'\ 



2 cos2 9 



\S\+0{R). 



The previous paragraph shows that there is an periodic 5'- 
DCC of density 5 = Ij^/n, where n = (2/\/3) + o(s^). So 
Theorem [3T| now implies that there exists a distinct difference 
configuration in S' containing at least m dots, where 



^(./2-2e- 



sin26') 



-R-o{R). 



V2VT- 



cos 9 

As in Subsection lVI-CI we may take 9 w 0.41586 to maximise 
this expression. Hence we have proved the following theorem: 

Theorem 34. Let fi « 1.61589 be the constant defined above 
TheoremU^ There exists a DD* (m, r) in which the number of 
o{R) « 0.86819r. 



dots is at least J -^^R 
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Fig. 12. Rotated square intersecting a circle 

TABLE I 

Upper and lower bounds on the number of dots in a distinct 
difference configuration 





lower bound 


upper bound 


[18] 


DD{m,r) 


{1/V2)r - o(r) 


(1/V2)r + 0(r2/3) 




DD(m,r) 


0.80795r - o(r) 


0.88623r + 0(r2/3) 


[19] 


T5T5*{m,r) 


0.79444r - o(r) 


0.86603r + 0(r2/3) 




DD*(m,r) 


0.86819r - o(r) 


0.95231r + 0(r2/3) 


[20] 



VII. Conclusion 

We introduced the concept of a distinct difference config- 
uration and gave specific examples for both the square and 
hexagonal grids for small parameters. We went on to provide 
general constructions for such configurations, as well as upper 
and lower bounds on the maximum number of dots such 
configurations may contain. In the case of distinct difference 
configurations using ]VIanhattan distance these bounds are 
tight asymptotically, as we have provided a construction for 
configurations which meets the leading term in our upper 
bound. For the remaining classes of configurations, there is 
a gap between the upper and lower bounds we have provided 
(see Table IJi. We believe the upper bounds to be realistic, and 
it is an interesting challenge to provide constructions that meet 
these bounds. 
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